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PRIORITY CLAIM 
This application claims priority of U.S. Provisional Application 
Serial No. 60/413,328, filed September 25, 2002. 

FIELD OF THE INVENTION 
The present invention relates generally to the field of hydraulic 
system analysis. 

BACKGROUND OF THE INVENTION 
Real-time health assessment for hydraulic pumps is a desired 
function due to, among other things, the high cost of unexpected failure of 
hydraulic systems. Typical hydraulic systems include both hydraulic- 
mechanical and electronic components, but most faults occur in the hydraulic- 
mechanical components. Common hydraulic system faults include, but are not 
limited to distortion, stress rupture, erosion, rubbing abrasion, impacting 
rupture, heat stress, and hot distortion. Furthermore, a hydraulic transmission 
and control system has its own special faults, such as oil pollution, leakage, air 
erosion, hydraulic blocking, pipe resonance, distortion of electrical signal, 
noise, and system surging. 

Many existing fault diagnosis methods for hydraulic systems are 
based on mechanical system parameters, with feature signals such as vibration, 
acoustic noise, and temperatures. However, because these parameters are 
indirect measures of hydraulic system operating conditions, and due to the 
multiple motion forms of hydraulic-mechanical components and the 
interference of multiple components via the hydraulic lines, a more 
complicated process is required to use these indirect parameters to monitor a 
state of the hydraulic system sensitively and accurately. 

For example, operation status of a hydraulic pump, a core 
component in a hydraulic system, directly influences the reliability of the 
hydraulic system. In piston-type hydraulic pumps, for example, common faults 



include, but are not limited to, worn pistons, swash plates, and distributing 
discs, bearing and spring failures, and loose piston shoes. These faults are 
often reflected in the pump discharge pressure, but are normally buried in the 
pulsation pressure signals. In addition, there are other noise sources, such as 
5 air erosion, hydraulic blocking, pipe resonance and leakage, etc. reflected in the 
pump discharge pressure. These noises normally result in a very low signal-to- 
noise ratio in the pump discharge pressure signals. Conventional health 
diagnosis methods, such as limit checking, spectrum analysis, and logic 
reasoning, require a distinguishable feature signal to detect faults, but these 
10 methods heretofore have not been sensitive or robust enough to reliably detect 
pump faults. 

To obtain more reliable pump health diagnosis results, vibration 
analysis methods based on spectral analysis have been disclosed. In an 
exemplary vibration-based diagnosis method, an accelerometer is installed on 

15 the shell of the pump, and fault diagnosis is performed by spectral analysis of 
the shell vibration signals. Diagnosis methods may include, for example: (1) 
calculating an over-limit mean square amplitude of the vibration, in which a 
fault state is diagnosed in the mean square value exceeds a preset threshold; (2) 
characteristic frequency analysis, in which the frequency spectrum of obtained 

20 vibration signals is compared with that of a normal vibration signal, where the 
fault signal and/or pattern is identified based on differences between the 
evaluating spectrum and the normal spectrum; and (3) time- frequency domain 
analysis, in which feature patterns are extracted based on signal distributions on 
both time and frequency domain signals to diagnose faults of the system. 

25 

SUMMARY OF THE INVENTION 
The present invention provides a method and apparatus for 
analyzing a hydraulic pump in real-time. A pressure signal is provided 
representing a discharge pressure of the hydraulic pump, and the pressure 
30 signal is decomposed into a plurality of levels. Each of the plurality of levels 
has at least one frequency band. A feature pressure signal is located in at least 
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one of the frequency bands and compared to a reference wavelet to determine if 
a fault exists in the hydraulic pump and/or a type of defect in the hydraulic 
pump. 

5 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a schematic of a hydraulic system having an apparatus 
for analyzing a hydraulic pump according to an embodiment of the present 
invention; 

10 FIG. 2 is a frequency-domain pulsating model of an axial piston 

pump; 

FIG. 3 is an illustration of a three-level wavelet decomposition of 
an original signal according to an embodiment of the present invention; 

FIGs. 4A - 4B together are a flowchart showing steps in a 
15 method for analyzing a hydraulic pump according to an exemplary embodiment 
of the present invention; 

FIGs. 5A — 5C together are a flowchart showing steps in a 
learning process for determining reference wavelets for diagnosis according to 
an exemplary embodiment of the present invention; 
20 FIG. 6 shows an experimental setup for testing a normal 

hydraulic pump and a defective hydraulic pump according to an embodiment of 
the present invention; 

FIGs. 7A-7D show experimental results for a measured pump 
discharge pressure from a normal pump, including an original pressure signal 
25 and wavelet coefficients for layers 1, 2, and 3, respectively, according to an 
exemplary method of the present invention; 

FIGs. 8A-8D show experimental results for a measured pump 
discharge pressure from a defective pump with loose piston shoes, including an 
original pressure signal and wavelet coefficients for layers 1, 2, and 3, 
30 respectively; and 
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FIGs. 9A-9D show experimental results for a measured pump 
discharge pressure from a defective pump with a worn swashing plate, 
including an original pressure signal and wavelet coefficients for layers 1, 2, 
and 3, respectively. 

5 

DETAILED DESCRIPTION OF THE INVENTION 

Fault diagnosis methods based on the vibration signals measured 

10 from the pump shell result in significant vibration effects, including 
environmental effects, which influence the quality of the obtained signals. 
More particularly, such signals consist of a series of harmonic frequencies and 
contain high background noises, making it difficult to distinguish feature 
signals of pump faults from the vibration signals. Furthermore, these methods 

15 are normally based on spectrum analysis that uses Fourier transform (FT), 
short-time Fourier transform (STFT), and/or time sequence analysis (TSA). 
Because these methods process signals in the frequency domain or time domain 
signals alone, they have limitations in practical applications. 

The present invention provides, among other things, a fault 

20 diagnosis method and apparatus that assesses, in real-time, pump health 
conditions and fault symptoms based on pump discharge pressure. 
Accordingly, preferred embodiments of the present invention can accurately 
predict a possibility of pump failure before such failure occurs, thus 
substantially reducing the likelihood of unanticipated hydraulic equipment 

25 failure and resulting downtime. By diagnosing and correcting a fault before it 
worsens to the point of pump and/or system failure, more expensive repairs or 
maintenance may be reduced. A lifetime of a pump and associated hydraulic 
system can be predicted based on a diagnosed fault. Reliability of systems 
having hydraulic pumps can be improved. 

30 A preferred method of the present invention analyzes a hydraulic 

pump in real-time by measuring discharge pressure of the pump, and applies 
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wavelet analysis to the measured discharge pressure to diagnose a fault. 
Generally, the wavelet analysis decomposes a pressure signal into a number of 
windows for evaluation, and compares one or more feature pressure signals 
within those windows to one or more reference wavelets. Reference wavelets 
5 are selected standard wavelets that correspond to a normally functioning 
hydraulic pump and/or hydraulic pumps having specific faults. By comparing 
the feature pressure signals to the reference wavelets, determinations can be 
made regarding the condition of the pump and the hydraulic system. 

More particularly, the preferred method provides a pressure 

10 signal representing the discharge pressure of the pump, and decomposes the 
pressure signal into a number of levels. Each of the levels has at least one 
frequency band. At least one feature pressure signal within one of the 
frequency bands is located and compared to at least one reference wavelet. The 
reference wavelet relates to a certain characteristic, and comparing the feature 

15 pressure signal and the reference wavelet in particular frequency bands can 
determine whether the characteristic exists in the hydraulic pump. By directly 
measuring the pump discharge pressure, environmental noise is significantly 
reduced. 

A preferred embodiment for analyzing a hydraulic pump includes 
20 a pressure sensor in fluid communication with a discharge port of the pump. A 
processor is coupled to the pressure sensor, and is configured to decompose a 
pressure signal from the pressure sensor into a number of levels, where each of 
the levels has at least one frequency band. The processor is configured to 
compare one or more feature pressure signals within at least one of the 
25 frequency bands to one or more reference wavelets to determine whether a 
characteristic exists. An alarm signal generator connected to the processor 
preferably is provided to generate an alarm signal indicating a pump fault, if a 
pump fault is detected. This alarm signal may, for example, alert a user of 
potential pump problems, request an on-site inspection, or order a replacement 
30 pump or pump component. The processor may be associated with one or more 
computers for on-site and/or remote monitoring, processing, and/or analysis. 
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Preferably, an embodiment of the present invention can be 
integrated into existing hydraulic pumps or systems without significant 
hardware modification. A hydraulic system is also provided having a hydraulic 
pump and an apparatus for analyzing the hydraulic pump. 

Referring now to the drawings, FIG. 1 shows a hydraulic system 
200 having a hydraulic pump 210 and an analyzing apparatus 220 for the 
hydraulic pump according to an exemplary embodiment of the present 
invention. An exemplary type of the hydraulic pump 210 is an axial piston 
fixed-displacement hydraulic pump. However, many other types of hydraulic 
pumps may be used. In addition to the hydraulic pump 210, the hydraulic 
system 200 includes various other hydraulic, hydraulic-mechanical, and/or 
electronic components, such as, but not limited to, a number of fluid passages 
222, a check valve 224, a motor 226 for the hydraulic pump, and a load for the 
hydraulic system, such as a relief valve 228. 

The analyzing apparatus 220 includes a pressure sensor 230 
mounted onto or otherwise integrated into the hydraulic system 200 and placed 
in fluid communication with a discharge port 240 of the hydraulic pump 210 
(FIG. 2). The pressure sensor 230 directly monitors the discharge pressure 
from the hydraulic pump 210. For example, the pressure sensor 230 may be 
placed downstream of the check valve 224. For example only, the pressure 
sensor 230 may be an Omega PX01C1-200G5T sensor, though others may be 
used. The pressure sensor 230 may have a range, for example, of 0 - 10 MPa 
and a bandwidth of 0 - 20 kHz, though other ranges and bandwidths are 
possible. 

In a preferred embodiment, a signal from the pressure sensor 230 
is transmitted to a suitable processor 250 via a suitable communication path 
260, which may be wired or wireless, and may or may not be part of a larger 
network. The processor 250, which is configured to analyze the signal, may be 
embodied in a computer having a Peripheral Connection Interface (PCI) card 
254 for connecting to the pressure sensor 230. It is also contemplated that 
other types of connections, boards, or cards may be used, or alternatively that 



the processor 250 may be a stand-alone device configured to perform analysis 
of the provided signal according to the present invention. Preferably, the 
processor 250 samples the signal from the pressure sensor 230 at discrete 
times, as a non-limiting example 500 Hz, to provide a pressure signal 
representative of the discharge pressure of the hydraulic pump 210. 

According to a preferred method, analysis of the hydraulic pump 
210 can be performed by analyzing only the measured discharge pressure of the 
hydraulic pump. A reason why this is possible is explained with reference to 
an exemplary, non-limiting model of the hydraulic pump 210 in FIG. 2, shown 
by example as a piston hydraulic pump connected to the load 228. In FIG. 2, 
all of the parameters are presented in the frequency domain. The hydraulic 
pump 210 in this embodiment is a positive displacement pump, which fills a 
cylinder 250 with hydraulic fluid when a piston 252 retracts and discharges the 
pressurized fluid when the piston extends. The most common defects in the 
hydraulic pump 210 are worn swash plates, pistons/cylinders, and piston shoes, 
loose piston shoes, and spring and bearing failure. Other defects are possible 
as well. 

The defects of the hydraulic pump 210 are reflected within 
certain frequency bands of pulsation pressure of the discharged fluid. Because 
the pulsation of discharge pressure is closely related to the flow pulsation, the 
following discharge flow pulsation model serves as a base for an exemplary 
pressure pulsation analysis. 

= q s +A\smayt\- qi 

where, q ° is the pump discharge flow, q *\s the rational discharge 

flow (discharge flow in normal operation), q * is the average pump discharge 

flow, A is the amplitude of flow pulsation, and qi is total leakage of the 
hydraulic pump 210. 
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The pump leakage, a function of the pump discharge pressure, 
plays an important role in the dynamic behavior of the pump discharge flow 
pulsation, and is defined as follows. 



o 



^9,+-^-{Po-Po) 

2 Po (2) 

5 where q ' is the average pump leakage, Po is the pump discharge 

pressure, and Po is the average pump discharge pressure. Equations (1) and (2) 
result in the following discharge flow pulsation model. 

9o = Is - <ll 

= q s +A\smat\-q l -^{p a -p a ) 

2 Po (3) 

To describe flow variations about the mean, Equation (3) can be 

10 rewritten as: 



Mo =(<lo -<lo) 

= (q s -?,)-(?/ -?/) 

= Aq s -^Ap 0 

2P 0 (4) 

The resulting pulsation model indicates that the discharge 
pressure fluctuation is affected by the pump discharge flow pulsation and flow 
15 fluctuation frequency, as well as pump leakage. Rewriting Equation (4) results 
in the following equation. 



q > (5) 



Equation (5) indicates that the actual discharge flow rate q ° from 
the hydraulic pump 210 is a function of the pump discharge pressure p 0 , 

20 hydraulic damping z * , and discharge flow rate under rated condition q s . 
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Observations from manual pump health diagnosis found that a 
loose or damaged piston shoe would result in a drop in the actual discharging 
flow, and a worn or damaged distributing disc would result in increased 
internal leakage and lead to a change in the pump hydraulic damping. Because 
it is often difficult to measure q * and z * directly, the pump discharge pressure 
Po and the pump discharge flow rate q ° can be selected as an indirect 
measurement reflecting the changes in q * and z * (Eqn. (5)). Though a model 
of a specific type of hydraulic pump 210 has been described, the present 
invention is not to be limited to analyzing only this specific hydraulic pump 
type. 

Given the pressure signal representative of the discharge pressure 
of the hydraulic pump 210, a wavelet analysis method is used for fault 
diagnosis. Generally, a wavelet analysis method decomposes spectral signals 
such as the pressure signal into windows in different frequency bands, and uses 
reference wavelets to characterize feature pressure signals within these 
windows. By this approach, fault diagnosis is performed according to the 
present invention by determining one or more appropriate reference wavelets to 
extract features from the pressure signals. 

A general wavelet can be defined using the following equation: 
( \ 1 (t-b\ 

v b AO=-rv — (6) 

*sla \ a J 

and a continuous wavelet transformation is defined as: 

WT{b,a) = ^£y^±y{t)d t (7 ) 

where, a and b are the scale and shift factors of the wavelet 

function. 

The scale parameter a scales the dimension of the window and 
the shift parameter b shifts the signal transformation window. By increasing a, 
the wavelet function reduces the time window, and vice versa. Therefore, 
wavelet analysis is capable of adapting the window dimension according to the 
signal frequency band. Parameter b, on the other hand, indicates the location of 
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the wavelet window along the time axis. By adjusting both parameters a and b, 
an appropriate size and location time window results for accurate and 
consistent identification of characteristic fault signals. Such features of 
wavelet transform analysis can improve the sensitivity and robustness of 
spectral signal analysis based pump health diagnosis. 

A general procedure in performing wavelet analysis is to first 
select one or more reference wavelets, and then compare located feature 
pressure signals with the reference wavelets using translated and dilated 
versions of the located feature pressure signals via shifting and scaling. There 
are many kinds of wavelets with different features, such as, but not limited to, 
the Harr wavelet, the Daubechies wavelet, the Morlet wavelet, and others, that 
can be selected as reference wavelets. 

In fault diagnosis of the hydraulic pump 210, the wavelet 
transform is used to identify singularities within the original pressure signals. 
Normally, a spectral signal such as the pressure signal may contain both non- 
continuous and non-continuous differential singularities. The non-continuous 
singular signals have abrupt changes at some characteristic points, which result 
in signal discontinuities. The use of wavelet transform can easily locate such a 
singularity. Non-continuous differential singular signals have abrupt changes 
in the first-order differential of the original pressure signals. In such a case, 
wavelet decomposition is applied on sampled pressure signals to locate the 
singularity within certain frequency bands. Windows within these bands are 
used for evaluating the pressure signal. 

A fault will result in specific singularities within a certain band. 
These faults are shown as variances in a wavelet coefficient, which is a 
coefficient indicating a difference or similarity between a feature pressure 
signal within a frequency band and the reference wavelet. The specific 
singularities cause corresponding wavelet coefficients to exceed their modular 
limits. If the reference wavelet represents a normally functioning pump, for 
example, a wavelet coefficient for the hydraulic pump 210 having a fault will 
exceed the maximum amplitude for healthy equipment for at least one 
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frequency band. Hence, the fault can be detected and located via wavelet 
analysis. 

In a preferred embodiment, a discrete wavelet transform (DWT) 
is used, and the signal from the pressure sensor 230 is digitally sampled by the 
processor 250 at discrete time steps to provide the pressure signal. By the 
DWT approach, the parameters, a,b, in Eqn. (6) are replaced using discrete 
values: a = 2 m ,b = 2 m n, and the continuous wavelets, y/ Qyb {t), are replaced by 
some orthonormal discrete wavelets: 2~ ml2 y/{2~ m t-n) , where m,n are scale 
factors within an integral space Z. The signal, /(/), is represented using the 
sum of its smooth approximation (low-pass) and its detailed description (band- 
pass): 



WIq CO 



/(0= Z</>^, n (0>^ 0 , n (0+ S E</>vv n (0>vv n (0 

w=-°° m=-co n=-co (8) 

m--oo 

where, <p(t) is a scaling function, P mo f(t) is the coarser 
approximation of /(/) in scale m 0 , and D m f(t) represents the differences 
among dilations. 

When a signal satisfies the relationship 
^-1/(0 = ^/(0 + ^/(0, it implies that the signal can be fine-scaled at 
P mo f(0 = fo and be decomposed into f 0 = P mo+l f(t)+D mo+l f(t) = f l +d x , where /, 
is the next coarser approximation of f 0 . The discrete model of wavelet 
analysis can therefore be represented as follows: 

fo = £<f, <P0,n (0 > <Po,n (0 = Z a 0,n<P0,n (0 
n--<X3 n=-oo 

/. = Z</.^(0>^(0= i«M^,..(0 (9) 

ft— —CO H--CO 

d t = Z< />,,,, (0>V,,«(') = £cd Un y Un «) 

n=-co n=-oo 

Using the same approach, f. can be further decomposed into 
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Based on this scheme, a set of examining signals such as pressure 
signals is decomposed using a low pass filter and a high pass filter, which 
results in two sets of sub-band signals. The sub-band signals are then 
reassembled to perform wavelet analysis. For example, when applying a three- 
5 level decomposition wavelet analysis to reassemble the original pressure signal, 
it will result in a wavelet coefficient vector, S, containing wavelet coefficients, 
a { and cdi , in both low and high frequency windows of these decomposed 
levels. A diagram illustrating a three-level decomposition in this manner is 
shown in FIG. 3, including, from top to bottom underneath the original signal, 
10 levels 1, 2, and 3, respectively. In harmonic analysis, such a decomposition 
procedure is called a 'two-channel 5 sub-band filtering scheme. The wavelet 
coefficient vector S can be displayed using the following expression. 

S = {a l9 cd x } 

= {{a 29 cd 2 } f cd t } (10) 
= {{{a 3 ,crf 3 },crf 2 },c£/ I } 

Preferably, to improve the analysis efficiency, the entire 

15 bandwidth of the pressure signal is normalized to be 1, so that, in the example 

of a three-level decomposition, the corresponding frequency band windows for 

evaluating the low frequency wavelet coefficients a, , a 2 and a 3 are 0-0.5, 

0-0.25 and 0-0.125, and the corresponding frequency bands for evaluating the 

high frequency wavelet coefficients cd l9 cd 2 and cd 3 are 0.5-1, 0.25-0.5 and 

20 0.125-0.25, respectively. As another example, for a ten-level decomposition, 

the first five high-frequency bands for determining high-frequency coefficients, 

would include .5-1, .25-.5, .125-.25, .0625-.125, and .03125-.0625, respectively. 

FIGs. 4A and 4B illustrate an exemplary, non-limiting method of 

diagnosing a fault in the hydraulic pump 210 using the pressure signal 

25 representative of the discharge pressure and DWT, according to the above 

description. The pressure signal is received by the processor (step 300), and 

may be tested, for example, by taking a covariance of the pressure signal, to 

determine if the pressure signal is reasonable. The pressure signal is 

decomposed into a number of decomposition levels (step 302). To decompose 
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the pressure signal, for example, the entire frequency band is set to be band ao 
at a level zero (step 304). The level "i" is set to zero (step 306), and then 
incremented (step 308). For each level "i" a low band filter, such as a digital 
filter embodied in the processor 250, filters the pressure signal (step 310) to 
produce a low-frequency band ai(step 311), and a high band filter, which also 
may be a digital filter, filters the pressure signal to produce a high-frequency 
band dj (step 312). When the desired number of decomposition levels is 
obtained (step 314), the processor 250 evaluates one or more of the frequency 
bands within windows using reference wavelets (step 316). It is possible, 
however, for evaluations of individual frequency bands to take place before the 
complete decomposition is completed. 

The number of decomposition levels needed for evaluation and 
the reference wavelets are determined according to a learning process, an 
example of which is illustrated in FIGs. 5A-5C. In the learning process, a 
default hydraulic pump is run having a known condition (such as a pump 
without defects, or a pump having particular, known faults) relating to a 
characteristic that is being evaluated in the tested hydraulic pump 210, and the 
pressure sensor 230 with the processor 250 provides a characteristic pressure 
signal representing the discharge pressure of the default hydraulic pump (step 
400). The characteristic pressure signal is decomposed (step 402) into a 
number of evaluation levels (decomposition levels used during the learning 
process, which may be more or less than the number of levels used during 
diagnosis) and corresponding frequency bands, in a similar process to the 
diagnosis process shown in FIG. 4A. The number of evaluation levels can vary 
according to processing time, etc., but it has been found that ten evaluation 
levels are typically sufficient. Less than ten levels may also be considered. 

In an exemplary, non-limiting method, the entire frequency band 
is set to be band ao (step 404). An evaluation level "i" is set to zero (step 406), 
and then incremented (step 408). For each evaluation level "i" the low band 
filter filters the pressure signal to produce a low-frequency band ^ (step 410), 



14 



and a high band filter, which also may be a digital filter, filters the pressure 
signal to produce a high-frequency band dj (step 412). 

After (or during) the decomposition, reference wavelets are 
selected, and feature pressure signals are identified as being similar to the 
5 reference wavelets. Individual feature pressure signals correspond to data point 
ranges within the decomposed frequency bands. Preferably, at least one 
reference wavelet is determined for each level "i" of decomposition. In an 
exemplary method, for one or more frequency bands of one or more 
decomposition levels, a feature pressure signal of the decomposed pressure 

10 signal within a particular data point range is identified that is similar to a 
standard wavelet (such as a particular Haar wavelet, Daubechies wavelet, 
Morlet wavelet, etc.). The standard wavelet chosen becomes the reference 
wavelet for that frequency band and, if only one frequency band is considered 
in a level, the reference wavelet for that level. 

15 In an exemplary, non-limiting method of identifying the reference 

wavelet as shown in FIGs. 5B - 5C, the level i is reset to zero (step 420) and 
incremented by one (step 422). Within each level, a number X, representing a 
particular possible reference wavelet (e.g., a particular standard wavelet), is 
reset (step 424) and incremented by one (step 426) to test each reference 

20 wavelet against candidate feature signals within a particular frequency band. 
The candidate feature signals are determined over a set of n data points within 
frequency bands. 

The candidate feature signal for a data point range n is compared 
to possible reference wavelet X (step 434) to determine a wavelet coefficient, 

25 which represents the difference between them. In some cases, the candidate 
feature signal has a similar proportional pattern to the possible reference 
wavelet, for example, but different amplitude. To provide an accurate 
comparison, since the pattern of the wavelet is the most significant detection 
tool, a candidate feature signal may be scaled (step 432) before comparing. In 

30 a non-limiting example, if a possible reference wavelet varies between 10 and - 
10 (a distance of 20), and a candidate feature signal varies between 2 and -2 ( a 



15 



distance of 4), each of the set of data points of the candidate signal wavelet is 
multiplied by a scaling factor of 5 for comparison with the possible reference 
wavelet. 

To compare the scaled candidate feature signal and the possible 
reference wavelet (step 434), similarities between candidate feature signals and 
possible reference wavelets are determined based on the wavelet coefficient. If 
the wavelet coefficient is substantially consistently within a relatively small 
band (for example, between -0.2 ~ 0.2), then the possible reference wavelet is 
selected as the reference wavelet. Since different possible reference wavelets 
are compared, the possible reference wavelet having the smallest wavelet 
coefficient band preferably is chosen as the reference wavelet. 

When all possible reference wavelets have been considered (step 
438), the reference wavelet and identified feature signal (i.e., a particular data 
range) are determined for the particular level (step 440). Once all levels have 
been considered (step 442), the learning process is completed. 

In the exemplary method in FIGs. 5A-5C, only candidate feature 
signals within the high frequency band di are compared to possible reference 
wavelets. However, it is contemplated that candidate feature signals within the 
low frequency band a; for each level may additionally or alternatively be 
compared and used as identified feature signals. In this case, there may be 
more than one reference wavelet for a particular level. 

Referring again to FIG. 4B, the frequency bands for the 
diagnosing process (step 316) preferably are evaluated at each level. In the 
non-limiting method of FIG. 4B, the level number i is reset (step 320) and 
incremented (step 322). For each level (and band, if more than one is evaluated 
in a level), wavelet transform is conducted on the pressure signal based on the 
reference wavelet for that level (and possibly band), to determine a wavelet 
coefficient. Put another way, the pressure signal data at a particular 
decomposition level (and possibly band) is converted to a wavelet coefficient. 

For example, the feature pressure signal (the pressure signal 
within an examining window having the same number of data points as the 
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reference wavelet) within a particular frequency band (as shown by example, 
the high frequency band di) of the decomposed pressure signal is located (step 
324) and compared to the reference wavelet for that level (step 325). The 
particular reference wavelet is determined by the learning process shown by 
5 example only in FIGs. 5A-5C. Based on the reference wavelet, a series of 
identified data points within the window will be identified within the frequency 
band as a data set. 

The identified data set in the extracted feature pressure signal 
preferably is scaled based on the reference wavelet, and the scaled data set is 

10 used to perform a wavelet transform to determine the wavelet coefficients. For 
example, the wavelet coefficient cdj represents a similarity or difference for a 
high-frequency band at decomposition level i. Alternatively, a wavelet 
coefficient a; may represent a similarity or difference for a low- frequency band 
at decomposition level i. 

15 Preferably, the wavelet coefficient is calculated so that the feature 

to be detected is present when the wavelet coefficient reaches or exceeds a 
certain threshold (step 304). A threshold can be established to determine 
whether a sufficient similarity or difference has been identified. For instance, 
if the wavelet coefficients of a normal signal are varying within a band c — c, 

20 the wavelet coefficients of a malfunction signal will exceed c. Therefore, c can 
be the threshold. If the wavelet coefficient meets or exceeds the threshold (step 
326), a determination of a fault is made, and an alarm signal may be produced 
(step 328) by a suitable alarm signal generator. By detecting the amount of 
similarity or difference, determinations can be made about the condition of the 

25 hydraulic pump 210. 

For example, if the reference wavelet represents a normal 
hydraulic pump without defects, than a wavelet coefficient that exceeds a 
threshold signifying a difference between the located feature pressure signal 
and the reference wavelet indicates that the hydraulic pump 210 is not 

30 operating normally, i.e. a fault exists. Furthermore, as particular faults have 
corresponding signature patterns at certain frequency bands, a located feature 
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pressure signal from a frequency band can be compared to a representative 
wavelet representative of that fault. 

In an example of detecting a particular defect, the learning 
process is repeated for a pump having a known defect, such as a worn swash 
5 plate. One or more reference wavelets are found for one or more levels. Thus, 
during the diagnosis process, a separate wavelet transform is performed on a 
decomposed pressure signal to determine the wavelet coefficients (for 
particular levels and/or bands) representing the known defect. A pump 
exhibiting the known defect may have, for example, a wavelet coefficient 

10 within a band of 0.2 ~ -0.2, where 0 equals complete similarity. If a threshold 
representing similarity between the located feature pressure signal of a pump to 
be tested and the reference wavelet exceeds a particular amount, for example, if 
the wavelet coefficient is within the band corresponding to the known defect, 
then a determination can be made that the hydraulic pump 210 exhibits the 

15 particular fault. 

If, as is preferred, wavelet coefficients are determined for a 
plurality of levels, evaluations can be made at one or more of the levels to 
detect whether a threshold exists. For example, three wavelet coefficients for 
three corresponding levels of decomposition can be detected. When the desired 

20 number of levels has been considered (step 330), the diagnosis process may 
repeat (step 332) as many times as desired to provide an ongoing real-time 
diagnosis. 

Accordingly, a preferred embodiment and method of the present 
invention can detect not only the presence of a fault in the hydraulic pump 210, 

25 but also the type or cause of the fault. In this way, appropriate action can be 
taken to prevent failure of the hydraulic pump 210 and/or the hydraulic system 
200 before it occurs. Preferably, the signal provider generates a signal (step 
328) if the threshold exceeds a particular value. By configuring the processor 
250, the provided signal can be analyzed for both the presence of a fault and 

30 the type of fault, if one is detected. 
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In an exemplary embodiment, original discharge pressure signals 
from a normal pump and two defective pumps were decomposed into high 
frequency windows of rf i , d * and d * using Haar wavelets as the reference 
wavelets. Diagnosis was conducted on a laboratory scale hydraulic pump 
health diagnosis research platform, as shown in FIG. 6. In a testing hydraulic 
system 600, two testing pumps, a normal pump 602 and a defective pump 604, 
were installed in parallel. The pressure sensor 230 was installed on the 
discharge port of each pump 602, 604 to collect the discharge pressure signals. 
An electrohydraulic servo control valve 610 was also installed, and another 
pressure sensor 230 was installed as well. When one of the pumps 602, 604 
was in testing, the other pump was shut off to avoid any possible inference to 
the discharge pressure of the testing pump. 

Under normal operating conditions, the pump discharge pressure 
will always have small fluctuations around its average pressure, and the 
variation of all wavelet coefficients within the high frequency bands considered 
should fall within normalized bands. FIG. 7 A shows the pump discharge 
pressure signal obtained from the normal pump 602, and FIGs. 7B-7D show the 
three-level high frequency wavelet coefficients cd u cd 2y cd 3 from the pressure 
signal (using Haar wavelets as the reference wavelets). These results showed 
that the discharge pressure from the normal pump 602 was stable and its 
pulsation amplitude was low. The wavelet analysis results indicated that the 
variation of all three wavelet coefficients were within the -1 to +1 range. 

FIGs. 8A-8D and 9A-9D show the test and analysis results 
obtained from a defective pump with loose piston shoes or a worn swash plate, 
respectively. Comparing the original signals shown in FIGs. 7A-7D, 8A-8D, 
and 9A-9D, the results indicated that the pulsations of the discharge pressures 
(FIGs. 7A, 8A, and 9A) from the pumps were very similar except for a slightly 
higher amplitude from the defective pumps. These results verified that the 
original pulsation pressure signals were not capable of providing sufficient 
information to support pump health diagnosis. 
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By comparison, the results from wavelet analysis (FIGs. 7B-7D, 
8B-8D, and 9B-9D) indicated that there were remarkable differences between 
the wavelet coefficients cd u cd 2 , cd 3 in all decomposed high frequency 
windows from defective pumps and from the normal pump. The resulting 
wavelet coefficients in all three windows from the normal pump exhibited 
relatively stable patterns, and the majority of the coefficient values were within 
the boundaries between -0.6 and +0.6 for cd * (FIG. 7B), between -0.4 and 
+0.4 for cd 2 (FIG. 7C), and between -0.3 and +0.3 for cd * (FIG. 7D). 

When a pump had loose piston shoes, the wavelet coefficients 
exhibited a harmonic pattern in all three layers. In addition, the amplitudes of 
these coefficients were also increased to between -0.8 and +0.8 for cd * (FIG. 
8B), between -0.7 and +0.7 for cd 2 (FIG. 8C), and between -0.6 and +0.6 for 
cd * (FIG. 8D). Such changes in the obtained higher wavelet coefficients 
clearly indicated a deviation from normal pump coefficients and therefore can 
be used to identify a pump defect. 

For the pump having a worn swashing plate, the wavelet 
coefficients did not show a harmonic pattern as had been seen from a pump 
with loose piston shoes. However, the amplitudes of these coefficients were 
consistently higher than those from the normal pump (between -0.9 and +0.9 
for cd * (FIG. 9B), between -0.8 and +0.8 for cd 2 (FIG. 9C), and between -0.5 
and +0.5 for cd > (FIG. 9D). The results obtained from both defective pumps 
indicated that the wavelet coefficients obtained from three high frequency 
windows changed when a different type of defect occurred. Furthermore, the 
patterns of the coefficient changes were different for different types of pump 
defects. 

As shown in this example, the original pulsation pressure signals 
(FIGs. 7A, 8A, and 9A) from the pumps 602, 604 were very similar, and thus 
were substantially uninformative for reliable health diagnosis for hydraulic 
pumps. However, by decomposing the provided signals into located feature 
pressure signals and comparing them with the reference wavelet, as shown by 
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the wavelet coefficients cd u cd 2 , cd 3 (FIGs. 7B-7D, 8B-8D, and 9B-9D), 
distinguishable changes can be found between wavelet coefficients for both the 
normal pump 602 and the defective hydraulic pump 604. These differences 
also provide distinguishable features that can be used to identify particular 
pump defects. Accordingly, the wavelet analysis method according to the 
present invention can improve the capability of diagnosing the health 
conditions of hydraulic pumps by decomposing the original pulsation pressure 
signals. Furthermore, the patterns and the amplitudes of wavelet coefficients 
obtained from different decomposed signal windows can be used to assess the 
types of hydraulic pump defects. 

While various embodiments of the present invention have been 
shown and described, it should be understood that other modifications, 
substitutions, and alternatives are apparent to one of ordinary skill in the art. 
Such modifications, substitutions, and alternatives can be made without 
departing from the spirit and scope of the invention, which should be 
determined from the appended claims. 

Various features of the invention are set forth in the appended 

claims. 
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